Signatures of spin-preserving symmetries in two-dimensional hole gases 
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We investigate ramifications of the persistent spin helix symmetry in two-dimensional hole gases in 
the conductance of disordered mesoscopic systems. To this end we extend previous models by going 
beyond the axial approximation for III-V semiconductors. We identify for heavy-hole subbands an 
exact spin-preserving symmetry analogous to the electronic case by analyzing the crossover from 
weak anti-localization to weak localization and spin transmission as a function of extrinsic spin-orbit 
interaction strength. 
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Control over spin relaxation is essential to the op- 
erational capabilities of spin-based semiconductor de- 
vices [l|,y|. A major advance in this respect has been the 
identification of an SU(2) symmetry that confines spin 
evolution to a characteristic topology denoted "persis- 
tent spin helix" (PSH) and prevents spin relaxation 0,3- 
The latter could be identified by means of optical ex- 
periments [J, |5[ in two-dimensional electronic systems 
with linear-in-momentum Bychkov-Rashba |6| and Drcs- 
selhaus [7[ type spin-orbit interactions (SOIs) (here de- 
pending on the parameters a and /3, respectively) of equal 
magnitude. As had already been suggested in Refs. [1,0, 
this symmetry at a = /3 also becomes manifest in the 
weak localization (WL) feature in magnetoconductance 
traces of disordered materials with finite SOI, as opposed 
to weak antilocalization (WAL) mediated by spin relax- 
ation [10| . Recent experiments confirmed theoretical pre- 
dictions that the WL signature persists for n-doped sys- 
tems even in the presence of nonnegligable intrinsic SOI 
that scales as k^ Ill4l3j. 



The question naturally arises whether the phenomenon 
of suppressed spin relaxation also exists in p-doped con- 
ductors. Here we hence investigate the generalization of 
the PSH symmetry arguments in the context of struc- 
turally confined heavy-hole (HH) states in III-V semi- 
conductors forming a two-dimensional hole gas (2DIIG). 
In these materials the spin is subject to strong SOIs 
which typically enhances spin relaxation. This feature 
is mainly attributed to the carrier density dependence 
of the spin splitting, that has been investigated ana- 
lytically by means of diagrammatic perturbation theory 
within the spherical approximation for one or more sub- 
bands [1J|. Other works consider weak (anti-) localiza- 
tion in hole gases based on a semianalytical [15[ as well 
as a scmiclassical and numerical [16[ treatment of 4 x 4 
Luttinger-Kohn models [l7[. Here we focus on strong 
confinement described by an effective 2x2 model of 
the HH ground state. Our treatment is not restricted 
to the spherical or axial approximations, which signifi- 
cantly widens the range of observable phenomena com- 
pared to prior models. The low dimensionality allows for 



the identification of relevant symmetries that are used to 
deduce optimum parameter regimes for controlling spin 
relaxation. The structure of our model is given by the 
Hamiltonian 



H — HkinO'o + ri2DHG ' "", 



(1) 



where Hkin denotes the kinetic energy, ctq the identity 
matrix, cr the vector of Pauli matrices, and n2DHG the 
effective spin-orbit field coupling to the spin. In con- 
trast to the corresponding expression ri2DEG for elec- 
trons, where fc-linear terms are dominant [9[, to lead- 
ing order ri2DHG is characterized by a cubic momentum 
dependence. This is in agreement with existing 2DHG 
models and results from coupling of the HH to the light 
hole (LH) subbands [H, [H]. The result ([1]) is obtained 



via a perturbative expansion of the standard Luttingcr- 
Kohn Hamiltonian [17[ in the basis given by the subband 
edge states in growth direction along the [001] crystal 
axis [20|. Our model is adequate to describe systems in 
which the energy level closest to the HH ground state is 
the LH ground state of the quantum well. This applies 
to typical zinc-blende structure materials, as can be in- 
ferred from their material properties and calculated band 
structures given e.g. in Ref. [20|. In Eq. ([1]) cr represents 
the subspace spanned by the HH states of spin angular- 
momenta ±3/2 and the hole spin-orbit field is given by 

J^2DHG = Punk (2) 

+ Ad {-7fc|fe + S[klx + kill - ik^ky{kyX + k^y)]] 
+ Ar {5kl{kyX + fc^y) + 7[-fcj]A - kly + ik^kyk]] 

with the intrinsic Dresselhaus parameters 
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and the Luttinger parameters 7 = (73 + 72)/2 and 
(^ = (73 ~ 72)72 as in Ref. [2l[. We denote the material 



„8u8u 
'^41 



and 6!f '' consistent with Ref. m. 



•^41 



constants Ck, 

By going beyond the axiahy symmetric situation, S = 
0, the above expression is a generahzation to previous 
models [18|, [lOl, [2^- This allows for the description 
of a broader range of materials and to avoid neglect- 
ing anisotropics that are of importance, for instance, 
in the plasmon spectra of HH systems [23|. Vertical 
confinement is modeled by a potential well with per- 
pendicular wavenumber (fc^) that displays a splitting 
Ahl = 272?!^ (fc^) /mo between HH and LH-bands. For 
further analysis the terms proportional to the small pa- 
rameter Ck are neglected in Eqs. pi4|) . since for a nar- 
row confinement they are dominated by the terms pro- 
portional to 61^^"" (fcy ) in realistic materials, as shown in 
Table 6.3 in Ref. Q. 

Furthermore, the linear Dresselhaus term ([3]) effec- 
tively rcscalcs the axially symmetric part of the cubic 
Dresselhaus contribution. Equation ([2]) results from se- 
quential perturbative expansions up to third order in k 
and to first order with respect to the inverse splitting 
Ahl""'^ and to E^ imposed on the crystal. The identi- 
fication of enhanced spin relaxation times in this work 
is closely connected with broken axial symmetry, since 
here a conserved quantity related to the spin degree of 
freedom can only be constructed in the presence of terms 
with both two- and threefold rotational symmetry in the 
extrinsic and the intrinsic SOI. Although our findings 
suggest that obtaining an exact PSH symmetry may be 
limited by the parameters of realistic systems, an approx- 
imate symmetry in the leading-order Fourier components 
of n2DHG causes a weakly perturbed crossover from WAL 
to WL, similar to electronic systems with cubic intrinsic 
SOI 



11 



The effect of the spin symmetry on the magnetocon- 
ductance G can be analyzed by organizing the transmis- 
sion entering the Landauer-Biittiker framework 24l. 125 1. 



h 






Tn + r, 



OD, 



(6) 

according to their spin quantum numbers a, a' in terms 
of diagonal spin-preserving channels Td and a spin off- 
diagonal contribution Tod- Here, cr, a' = ±1 refer to 
an arbitrary basis defined in the ballistic leads of a two- 
terminal device representing our numerical model, while 
n, 771 are integers that define the transverse channel of the 
in- and outgoing states due to a hard-wall confinement 
of the leads. The lead wavcfunctions \(i>n.a) and \<^m.a') 
enter into the Fisher-Lee relation for the amplitudes 



ina.ma' OC \ Qi^^^^^^^T^ {KA Yl) (YI I Gr |y2) (y2 \4>ra,o') , 

where the integration is taken over the lead cross sec- 
tions [Hi. Gr = (Sf -H + 0+)"Ms the Green's func- 
tion of the scattering region at fixed Fermi energy E-p. 

Knap et al. [9| found in n-type systems particular re- 
lations between extrinsic and intrinsic SOI magnitude, 
for which the Coopcron becomes separable and a WL 
signal rather than WAL is observed. In terms of the 
structure provided by Eq. ([6]), Tod vanishes in this case 
and correspondingly spin scattering is absent even in 
transport in disordered systems. This is equivalent to 
the observation that the system displays an exact, dis- 
order independent symmetry [^ |3|, which allows for a 
decomposition within the corresponding constant eigen- 
basis llxa)} into n-er = Y.a=±\ Ea{^) \Xa) (x-t|- Hence, 
when taking the spin trace in Eq. ([6]) in the basis 
{IXct) = (l;Crexp[±z7r/4] )'''}, corresponding to the exis- 
tence of the conserved quantity I]± — Og^rkoy or, equiv- 
alently, fixed in-plane spin orientation along 9 ~ ±t:/A^ 

one finds that Tqd oc J^a^a' I(x<t|x<t')I^ = J^a^a' ^<y-<y' 
is supressed and Td decomposes into two independent 
channels which trivially display WL [9|. 

In the hole model ([2]) we find the analogue to the elec- 
tronic PSH symmetry if the system parameters fulfill 
Ar/Ad = ±1 and 7 = —(5, i.e., 73 = 0. In these two 
cases the direction of r22DHG is fixed independently of 
the momentum, more precisely by r22DHG oc [— fc|,(fca;± 
ky) ± 3kxky{kx ± ky) — fc^ T k^]ix ± y). We illustrate 
these cases in Fig. [1] where the effective spin-orbit field 
^2DHG is oriented along a fixed direction for both spin- 
split subbands [3^. By engineering the quantum well 
width and choosing the electrical gate field appropriately, 
a regime of enhanced spin lifetime should be accessible 
in hole gases. In realistic material systems it is possible 
to influence the effective values of 73 [35| . In particular 
at Ar/Ad = — 1 and 73 = 0, as investigated in Fig. [21 the 
symmetric situation is obtained. 

The symmetry condition becomes apparent in the mag- 
netoconductance of disordered 2DHG systems, as illus- 
trated in Fig. [21 There we show results of the numerically 
calculated disorder averaged transmission, Eq. ^ , for fi- 
nite cubic intrinsic SOI Ad as a function of the extrinsic 
SOI Ar, respectively the parameter 77 = Ar/Ad, setting 
73 = 0. Representative examples of the conductance cor- 
rection traces are shown in Fig. [Ifa), which display typ- 
ical WAL and WL lineshapes as a function of magnetic 
flux (f> from a homogeneous magnetic field perpendicu- 
lar to the 2DHG plane. Considering the dependence on 
V — Ar/Ad, wc find pronounced signatures of WAL if 
77 is far from —1. When 77 approaches the value of —1, 
a crossover from WAL to WL occurs as indicated by a 
maximum negative conductance correction, in agreement 
with the symmetry argument. In Fig.[2ljb) our results are 
summarized in terms of the conductance at maximum 
magnetic flux (T((/)i„ax)) subtracted from the correction 
at zero flux {T{0)) plotted as a function of rj, where we 
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Figure 1: (Color online) Fermi surface for the different spin di- 
rections obtained from Eq. IT} (violet and blue contours) and 
the corresponding direction of the effective spin-orbit field, 
fi2DHG, illustrated by arrows. The SOI parameters establish 
a persistent spin helix for holes with uniaxial spin orientation 
corresponding to Ar/Ad — -1-1 (a) or Ar/Ad = —1 (b). In 
both cases the Luttinger-Kohn parameters are 7 = —(5, i.e., 
73 = 0. 




Figure 2: (Color online) Signatures of spin-preserving symme- 
tries in weak localization of a two-dimensional hole gas. (a) 
Disorder-averaged magneto conductance correction l^T{(j))) — 
{T{(f)ma,x)) as a function of fiux (f) (in units of 00 = h/e; 
4'ni!!Lx/(t>o = 3.1) for spin-orbit coupling ratios 77 = Ar/Ad ~ 
-1.65, -0.017, -1.33, -0.67, -1 (from top to bottom), 
(b) Conductance correction (r(0)) — (r((^max)) as a function 
of 77. Negative magneto conductance reflects suppression of 
spin relaxation close to r; = — 1. System parameters used in 
(a) and (b): Disorder average over 1000 impurity conflgura- 
tions for a scattering region of aspect ratio (length to width 
W) 200:80 unit cells with periodic boundary conditions in 
transverse direction. Quantum transmission computed for 
tpW/TT = 13 hole states per spin supported in the leads, elas- 
tic mean free path I = 0.04iy, 73 — 0, and fixed Dresselhaus 
spin precession length knW ^ 1, defined below Eq. ((8]). 



chose ( 



3.100- The results show that the parameter 



regime where a PSH type symmetry occurs is character- 
ized by a negative conductance correction, i.e., by a WL 
signature. 

Our analysis above is confirmed within a diagram- 
matic perturbative treatment by exact diagonalization 



of the Cooperon C'(Q) in the framework of the effective 
model (fll) . For this purpose the scheme presented in 
Refs. [23, UM fo'^ electrons is generalized to holes. The 
diagrammatic approach is justified since we assume to be 
in the diffusive regime, EpT 3> 1, with elastic scattering 
time T. Since we consider the spin ±3/2 subspace, one is 

left with ii 



C(Q) 



Dh {hQ + ^m„ (a) ■ S) + H, 



(7) 



where D^ = tu|/2 is the diffusion constant, Q = k + k' 
the Cooperon momentum composed of momenta of the 
retarded and advanced propagators with spins 3/2?icr and 
3/2;icr', giving rise to S = 3/2h{(T + a'). In Eq. ^ the 
average over all directions of the velocity vp on the Fermi 
surface is denoted by (...), and a is defined via the re- 
lation (T • r22DHG ^ k ■ (a ■ (j) describing the Rashba and 
Dresselhaus SOI. H,. is a Q-independent term which gen- 
erally leads to spin relaxation. However it can be shown 
that one of the gaps of the eigenmodes in the triplet sec- 
tor of the Cooperon vanishes if the aforementioned sym- 
metries are present. Since the Cooperon gaps can be in- 
terpreted as spin-relaxation rates [27|, |28[ , this supports 
the numerical findings of persistent spin states. 

The numerical setups on which our simulations are 
based are disordered hole systems connected to two ter- 
minals, represented by ballistic semi- infinite leads with- 
out SOI. The latter is switched on and off adiabatically 
over one fifth of the total length of a rectangular scat- 
tering region to which the leads are connected. We use 
an average over an Anderson-like uniformly distributed 
random-box potential Vdis to simulate disorder. The per- 
pendicular magnetic field is included by means of the 
Peierls substitution. The Hamiltonian is then discretized 
on a tight-binding grid in position space and the trans- 
mission amplitudes are obtained by an optimized recur- 
sive Green's function algorithm [30|. Since we are inter- 
ested in modeling bulk transport, we implemented peri- 
odic boundary conditions in the transverse direction to 
minimize effects from the boundaries. 

Apart from considering the indirect influence of the 
PSH symmetry on the WL-WAL transition, it seems nat- 
ural to search for a manifestation of a symmetry in Td, 
Eq. ^ , since its effects could be determined by magnetic 
polarization of the leads, allowing for spin transistor op- 
eration even in the presence of disorder [2] . Numerically 
we can confirm the validity of the latter approach by cal- 
culating the normalized quantity Ty)/(Td + Tod) as a 
function 77 = Ar/Ad, as shown in Fig. [31 We identify 
a pronounced transmission maximum at 77 = 1 in the 
basis corresponding to the -|-7r/4 spin orientation even 
in situations where the exact PSH-type symmetry is not 
realized. In the given example we chose the Luttinger- 
Kohn parameters 72 = 1 and 73 = 0.25 which correspond 
to such a nonideal setup. For parameters far from 77 = 1 




Figure 3: Ratio of disorder averaged diagonal transmission 
over total transmission (Tu) / (T) = (Td> / (To + Tod) as a 
function of rj for a scattering region with 150:80 aspect ratio 
for fixed Dresselhaus spin precession length fc^^ = 1.3 W and 
Luttinger-Kohn parameters 72 = 1 and 73 = 0.25 for which 
no exact spin-preserving symmetry can be established. The 
peak of {To) / (T) at r; = 1, coincides with the maximum of 
the diabatic transition probability of Eq. ((8]) (indicated by the 
arrow). Average Transmission shown includes 1000 disorder 
configurations. 



the spin transmission is equally distributed among the 
diagonal and ofF-diagonal channels. 

When \ri\ approaches 1, Td formally corresponds 
to the probability of diabatic Landau-Zener tran- 
sitions between instantaneous eigenstates |±J^) = 
(1, ± exp[— i arctan {fly/^x)])^ of the spin-orbit contribu- 
tion ([2]). The momentum direction is changed by disor- 
der scattering such that the spin evolution is subject to 
inhomogeneties of the effective spin-orbit field fi. 

At the minima of the anisotropic spin splitting 
2\Ecr{Cl)\, this induces transitions of the type |±ri) — > 



31 



J^ri) with Landau-Zener transition probability Pd 
|32| • These transitions enhance the value of Td while com- 
pletely suppressing Tqd for Pd = 1- The spinors {|xcr)} 
underlying Eq. (JG)) coincide with the diabatic superposi- 
tion of the states |±r2). The latter can be checked by 
considering {xa\^ ■ ""IXo-)- Within the HH model ([2]) 
the diabatic basis coincides with that of the PSH eigen- 
states IIxct)} of a 2DEG [3- For p-type systems we find 
a probability [36[ 



ln(Pi 



DJ2DHG 



Cl\kn\\i + S\{l-\7j\Y 



(8) 



with the elastic mean free path / and a correction factor 
C entering together with the transport time t into the 
rate of change in angle 9 in the relation 69 = TT/2St/{TC,). 
Equation ([5]) is derived under the assumption that 7 ^ 
—(5. Although the expression for the Landau Zener tran- 
sition probability predicts a clear maximum at |?/| = 1, 
Equation ([8]) does not cover the description of Td for pa- 
rameters where the PSH symmetry is established. It is 
nevertheless applicable to realistic material parameters 
when 73 7^ and consequently 7 7^ (5, which is verified 
by a numerical transport analysis. 

The analyis of Td can be applied to electronic systems 
as well, with an effective spin-orbit field, 

^2DEG = akxz + l3{kxX - kyy)+j {-k^-kyX + kykly) , 

(9) 



for transport along the [100] direction in a 2DEG grown 
in [GDI] direction and with (k'j) ^ = P M- I'^ systems 
described by this model the corresponding Landau-Zener 
transition probability is given by 



ln(PD)2DEG = Cl\kp\ (r^/2 -l±r,f 



(10) 



with the Dresselhaus spin precession length kg = 

(TOeff/3/fi'^) J ratio of cubic and linear SOI F^ = fc|7//3 
and correction factor ( as it appears in Eq. ([8|). This 
model has been verified by numerical calculations which 
are beyond the scope of this work. 

In both p- and n-type systems, the signatures in Td are 
robust against disorder. Therefore as an experimental 
approach to analyzing spin relaxation lengths in trans- 
port within HH systems, a detection of the PSH signa- 
ture in the longitudinal conductance of a spin-polarized 
current is favorable. The mechanism responsible for the 
peaks in Td the momentum space analogue to the ef- 
fect of a spatially inhomogeneous helix-type Zeeman term 
on the spin conductance in dilute magnetic semiconduc- 
tors [331 1 . An alternative measurement method for further 
investigation of the HH PSH is represented by magneto- 
optical Kerr rotation techniques, which recently allowed 
to map the spin topology in electronic systems [5[. 
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